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We discuss probability distributions for the cosmological constant A and 
the amplitude of primordial density fluctuations Q in models where they both 
are anthropic variables. With mild assumptions about the prior probabili- 
ties, the distribution -P(A, Q) factorizes into two independent distributions 
for the variables Q and y oc A/Q^. The distribution for y is largely model- 
independent and is in a good agreement with the observed value of y. The 
form of P{Q) depends on the origin of density perturbations. If the perturba- 
tions are due to quantum fluctuations of the inflaton, then P{Q) tends to have 
an exponential dependence on Q, due to the fact that in such models Q is 
correlated with the amount of inflationary expansion. For simple models with 
a power-law potential, P{Q) is peaked at very small values of Q, far smaller 
than the observed value of 10^^. This problem does not arise in curvaton-type 
models, where the inflationary expansion factor is not correlated with Q. 



I. INTRODUCTION 

The fundamental theory of nature may admit multiple vacua with different low-energy 
constants. This possibility has attracted much attention in the context of "string theory 
landscape" [1-3]. When combined with the theory of eternal inflation [4,5], it leads to 
the picture of a "multiverse" , where constants of nature take different values in different 
post-inflationary (thermalized) regions of spacetime. The key problem in this theoretical 
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framework is to calculate the probability distribution for the constants. Once we have the 
distribution, we can use the principle of mediocrity - the assumption that we are typical 
among the observers in the universe - and make predictions for the constants at a specified 
confidence level. 

A major success of this program has been the prediction of a non-zero cosmological 
constant A [6-9]. The observed value of A is well within the 2a range of the theoretical 
distribution. Recently, however, it has been argued that this successful prediction does not 
survive when other parameters, such as the amplitude of primordial density fluctuations 
Q, are also allowed to vary [10,11]. Here, we critically examine these claims and discuss 
some issues surrounding the calculation of prior probabilities. We start with a review of the 
original prediction for A. 

After this work was completed, we noticed the paper [12] by Feldstein, Hall and Watari, 
which has a substantial overlap with our work. 

II. ANTHROPIC PREDICTION FOR A 

The key observation, due to Weinberg, is that structure formation effectively stops when 
A comes to dominate the universe. If no galaxies are formed by that time, then in the future 
there will also be no galaxies, and therefore no observers. The corresponding anthropic 
bound on the cosmological constant is [6] 

A < 500p^o (1) 

where p^o is the present matter density. A lower bound on A can be obtained by requiring 
that the universe does not recoUapse before life had a chance to develop. Assuming that 
the timescale for life evolution is comparable to the present cosmic time, one flnds [13,14] 
A ^ — PmO- 

The anthropic bound (1) specifles the value of A which makes galaxy formation barely 
possible. However, if the effective A varies in space, then most of the galaxies will not be 
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in regions characterized by this marginal value, but rather in regions teeming with galaxies, 
where A dominates after a subatantial fraction of matter has already clustered. 

To make this quantitative, we define the probability distribution P{A)dA as being pro- 
portional to the number of observers in the universe who will measure A in the interval dA. 
This distribution can be represented as a product [7] 

P{A)dA^nobs{A)Pprior{A)dA. (2) 

Here, Pprior{A)dA is the prior distribution, which is proportional to the volume of those parts 
of the universe where A takes values in the interval dA, and nobs (A) is the number of observers 
that are going to evolve per unit volume. The distribution (2) gives the probability that a 
randomly selected observer is located in a region where the effective cosmological constant 
is in the interval dA. 

Of course, we have no idea how to calculate riobs, but what comes to the rescue is the 
fact that the value of A does not directly affect the physics and chemistry of life. As a rough 
approximation, we can then assume that nobs (A) is simply proportional to the fraction of 
matter / clustered in giant galaxies hke ours (with mass M ^ Mq — lO^^Mo), 

no6s(A) oc /g(A). (3) 

The idea is that there is a certain number of stars per unit mass in a galaxy and certain 
number of observers per star. The same approximation can be used for other "life-neutral" 
parameters, like the amplitude of primordial density fluctuations Q or the neutrino masses. 

Now we have to decide what should be used for the prior probability in (2). At present, 
the details of the fundamental theory and of the inflationary dynamics are too uncertain for 
a deflnitive calculation of the prior. We shall instead rely on the following general argument 
[15,16]. 

Suppose some parameter X varies in the range /S.X and is characterized by a prior 
distribution Pprior{X). Suppose further that X affects the number of observers in such a 
way that this number is non-negligible only in a very narrow range ^Xobs ^ AX. Then one 

3 



can expect that the function Pprior{X) with a large characteristic range of variation should 
be very nearly a constant in the tiny interval A.Xobs- In the case of A, the range AA is set 
by the Planck scale or by the supersymmetry breaking scale, and we have (AA)o6s/AA ~ 
10-60 _ 10-120 Hence, we expect 

Ppriori^) ~ const. (4) 

We emphasize that the assumption here is that the value A = is not in any way special, 
as far as the fundamental theory is concerned, and is, therefore, not a singular point of 

-^rior (A) . 

Combining Eqs.(2),(3),(4), we obtain 

P(A)(x/g(A). (5) 

The fraction of clustered matter, /g(A), can be calculated using the Press-Schechter approx- 
imation. Restricting attention to positive values of A, one finds [9,17] 



/g(A) oc erfc 



0.8 



A 



(6) 



where ac is the density contrast on the galactic scale and pm is the matter density. The 
product PmO'G time-independent during the matter era, so it can be evaluated at any time. 
Introducing a dimensionless variable 

y = A/p^4, (7) 

we can write the distribution as 

dP{y) oc erfc(0.8y^/^)yci(lny). (8) 
This distribution is peaked at y ~ 1. Discarding 2.5% at both tails of (8) yields the range 

0.043 <y<16. (9) 
This is the prediction for y at 95% confidence level. 



The product Pm'^G expressed in terms of the present-day values, p^o and ctgo- 

Hence, y is & measurable quantity, and the anthropic prediction can be tested. As already 
mentioned, the observed value of y is well within the 95% range (9) (see, e.g., [18]). 

Apart from this successful prediction, two additional factors make the anthropic expla- 
nation of the observed value of A particularly compelling. First, there are no plausible 
alternatives. Second, the anthropic approach also gives a natural resolution to the coinci- 
dence puzzle: Why do we happen to live at the very special epoch when A ~ If we 
denote the redshifts at the epochs of galaxy formation and of A-domination by zq and za, 
respectively, then most of the galaxies should be in regions where zj^ ^ zq- Regions with 
z\ ^ Zg will have very few galaxies, while regions with zp^ <^ zg will be rare simply because 
they correspond to a very narrow range of A near zero (small "phase space"). Hence, we 
expect a typical galaxy to be located in a region where 

za ~ Zg- (10) 

The galaxy formation epoch, ~ 1 — 3, is close to the present cosmic time. This explains 
the coincidence [19,20]. 

III. VARIABLE Q 

Banks et. al. [10] and Graesser et. al. [11] have argued that the successful prediction 
for A will be destroyed if the density fluctuation amplitude Q is also allowed to vary. Q is 
defined as the density contrast at the time of horizon crossing. It is approximately the same 
on all scales of astrophysical interest. The anthropically allowed range of Q is [21] 

10-^ < g < 10-^ (11) 

The observed value, Q ~ 10~^, is in the middle of this range. 

For larger values of Q, galaxies form earlier, so A can get larger without interfering 
with the galaxy formation process, h.jnax oc Q^. Since larger A correspond to larger phase 
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space, this suggests that the probabihty is maximized for large Q and A, e.g., Q ~ 10 ^, 
A ~ lOVmo- 

To examine the situation more carefully, let us now consider the probability distribution 
for A and Q. Since the galactic-scale density contrast (Tg is linearly related to Q, ug oc Q, 
we can equivalently consider the distribution for A and ctg, 

dP(A, ac) a Pprior(A, aG)eM^-^y^'^)d^d<^G, (12) 

where, as before, y is given by Eq. (7). The narrow anthropic range of A suggests, also as 
before, that the prior probability is independent of A in the range of interest, 

PprioriK (^g) Ppri(yr{(^ g) ■ (13) 

Substituting this in (12) and changing variables from (A, ctg) to {y,(JG), we obtain [19] 

dP{y, aa) oc a%Ppri<yr{(yG)d(TG ■ erfc(0.8y^/2)dy. (14) 

Remarkably, the distribution for y and ctg is factorized. Moreover, the y-part of the 
distribution is exactly the same as it was for the original model, where ug was not variable. 
This means that the successful prediction for y is completely unaffected by the variation of 
ac [17]. 

IV. THE LARGE Q CATASTROPHE 

Let us now consider the distribution for ctg, 

dP{aG) oc aGPprior{o-G)dcrG. (15) 

The value of aG is determined by the horizon-crossing perturbation amplitude Q, which is 
in turn determined by the inflaton potential V{(j)), 

y3/2 

(JgocQ (X ——, (16) 
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where the right-hand side is evaluated at the value of corresponding to horizon crossing 
for the comoving scale of interest. With a power-law potential 

v(<t>) = xr, (17) 

this gives, neglecting logarithmic factors (see, e.g., [22]), 

Q oc AV2. (18) 

Suppose now that A is a variable, which may be determined by some additional scalar 
field. It seems natural to assume that the range of variation for A is A A ~ 1. The ob- 
served value of aa is obtained for A ~ 10~^^, and the anthropic range (11) corresponds to 
10~^^ ^ A ^ lO"-*^^. Since this range is so narrow, Graesser et. al. [11] argue that the same 
logic we used for the cosmological constant implies that the prior for A should be nearly flat 
in the range of interest, 

dPprior/dX ^ const. (19) 
Then it follows from (18) that dPprior{Q) cx: QdQ, and Eq. (15) gives 

dP{Q) oc Q^dQ. (20) 

This distribution is strongly biased in favor of large values of Q. If anthropic factors cut 
off the distribution above ^ 10~^, then Eq. (20) suggests that this cutoff value is 10^ 

times more probable than the observed value Q ~ 10^^. This is the large Q catastrophe. 

We note, however, that there is an important difference between the cosmological con- 
stant and the inflaton self-coupling A, which may invalidate the argument for the flat prior 
(19). Unlike the small cosmological constant, the value of A has a strong effect on the dynam- 
ics of inflation. As we shall see later, smaller values of A give larger inflationary expansion 
factors. Hence, A = is a rather special value of the coupling, and the flat prior assumption 
is not justified. 
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V. THE PROBLEM OF THE PRIOR 



The calculation of prior probabilities is a very challenging and important problem. The 
main difficulty is that the volume of thermalized regions with any given values of the pa- 
rameters (such as Q or A) is infinite, even for a region of a finite comoving size. To compare 
such infinite volumes, one has to introduce some sort of a cutoff. For example, one could 
include only regions that thermalized prior to some time tc and evaluate volume ratios in 
the limit tc ^ oo. However, one finds that the results are highly sensitive to the choice of 
the cutoff procedure (in this case, to the choice of the time coordinate t ) [23]. (For a recent 
discussion, see [24,25].) 

An eternally inflating universe consists of isolated thermalized regions embedded into 
the inflating background of false vacuum. These thermalized islands are rapidly expanding 
into the inflating sea, but the gaps between them are also expanding, making room for more 
thermalized islands to form. The thermalization surfaces at the boundaries between inflating 
and thermalized spacetime domains are three-dimensional, infinite, spacelike hypersurfaces. 
The spacetime geometry of an individual thermalized domain is most naturally described by 
choosing the corresponding thermalization surface as the origin of time. The domain then 
appears as a self-contained, infinite universe, with the thermalization surface playing the role 
of the big bang. Following Alan Guth, we shall call such infinite domains "pocket universes" . 
All pocket universes are spacelike-separated and are, therefore, causally disconnected from 
one another.^ 

The reason for the cutoff-dependence of probabihties is that the volume of an eternally 
inflating universe is growing exponentially with time. The volumes of regions with all pos- 
sible values of the parameters are growing exponentially as well. At any time, a substantial 
part of the total thermahzed volume is "new" and thus close to the cutoff surface. It is not 



^In models where false vacuum decays through bubble nucleation, the role of pocket universes is 
played by individual bubbles. 
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surprising, therefore, that the result depends on how that surface is drawn. 

It has been suggested in [26] that the resolution of this difficulty may lie in the direction 
of switching from a global distribution, defined with the aid of some globally defined time co- 
ordinate, to a pocket-based distribution. In a particular case, when the variable parameters 
take all their values within each pocket universe, the problem appears to have been solved 
completely. In this case, all pocket universes are statistically equivalent, and we can pick 
any one of them to calculate the distribution. The prior probability for the constants can 
be identified with the volume fraction occupied by the corresponding regions of the pocket 
universe. We can find this fraction by first evaluating it within a sphere of large radius R 
and then taking the limit it! — > oo. This is the so-called "spherical cutoff method". 

In general, however, there are several distinct types of pocket universes, and we have to 
face the problem of comparing volumes in different pockets. It is not clear how this can 
be done, since the spherical cutoff method cannot be applied to disconnected spaces. The 
problem of defining probabilities in the general case still remains unresolved. A conjecture 
toward this goal will be discussed in the following section. 

To get some idea of what the prior distribution for Q may be like, we shall first focus on 
the case of identical pockets, where we know how to calculate probabilities. The spacetime 
structure of a pocket universe is illustrated in Fig. 1. The surface S^, in the figure is the 
thermalization surface. It is the boundary between inflating and thermalized domains of 
spacetime, which marks the end of inflation and plays the role of the big bang in the pocket 
universe. The surface is the boundary between the stochastic domain, where the dynamics 
of the inflaton field is dominated by quantum diffusion, and the deterministic domain, where 
it is dominated by the deterministic slow roll. Thus, marks the onset of the slow roll. 
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quantum diffusion 



FIG. 1. The spacetime structure of a pocket universe. 

The prior probabihty Pprior{X)dX for a variable X is defined in terms of the volume 
fraction on the thermalization hypersurface E*. It can be expressed as 



where Pq{X) is the distribution (volume fraction) on S^, and Z{X) is the volume expansion 
factor during the slow roll. We assume for simplicity that the field responsible for the value 
of X interacts very weakly with the infiaton, so that X does not change appreciably during 
the slow roll. Otherwise, we would have to include an additional Jacobian factor in Eq. (21), 
as indicated in [26]. 

The expansion factor Z{X) can be found from 



Pprior{^) OC Pq{X)Z^ 



in 



(21) 




(22) 



where 



H{cf>;X)^[8nV{cf>,X)/3] 



1/2 



(23) 
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is the inflationary expansion rate, V{(f), X) is the inflaton potential, and prime stands for 
derivative with respect to 0. and ^^^x are the values of (p at the boundary surfaces Eg 
and E*, respectively. They are defined by the conditions 

H'/H\<Pgx,X)^l (24) 

and 

H'/H{4>,x,X)r^l. (25) 

The distribution Pq{X) can in principle be determined from numerical simulations of 
the quantum diffusion regime. Some useful techniques for this type of simulation have 
been developed in [23,27]. The analysis is greatly simphfied in a class of models where the 
potential is independent of X in the diffusion regime, and Pq{X) can be determined from 
symmetry, 

Pq{X) = const. (26) 

An example is a "new" inflation type model with a complex inflaton, = |0|exp(iX). 
Inflation occurs near the maximum of the potential at = 0, and we assume in addition 
that the potential is symmetric near the top, V = V{\(j)\). Eq. (26) follows if this property 
holds throughout the diffusion regime. In such models, the distribution (21) reduces to 

Pprior{X)^Z\X). (27) 

Eq. (27) has a simple intuitive meaning: the prior probability is determined by the volume 
expansion factor during the slow roll. 

VI. A CONJECTURE FOR A MORE GENERAL CASE 

In more general models, the factor Pq{X) may provide an additional dependence on X. 
It has been recently conjectured [28] that this factor is of the form 
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Pg{X)(xeMS{X)], 



(28) 



where 



S{X)=7T/H'{cj>,x,X) 



(29) 



is the Gibbons-Hawking entropy of de Sitter space. 

The idea is that the evolution in the diffusion regime is ergodic, in the sense that all 
quantum states are explored with equal weight. The conjecture is that this ergodic property 
extends to the diffusion boundary Hg. The right-hand side of Eq.(28) is simply the number 
of quantum states for the specified values of the parameters (fields) X on the boundary. The 
conjecture (28) has been verified in some specific models, but the conditions of its validity 
are not presently clear. 

In models with pockets of several different types, we need to introduce an additional 
factor pj, characterizing the relative probability of different pockets. The full expression for 
the prior probability is then given by 



where the index j labels different types of pockets. A natural choice for the factor pj is the 
so-called comoving probability, which can be defined as the probability for a test particle, 
starting near the top of the potential V{(f),X), to end up in a pocket of type j. This 
probability can be easily calculated using the Fokker-Planck equation. 

One problem with this definition of pj is that it has some dependence on the initial 
conditions assumed for the test particle. For example, if the potential has several peaks of 
comparable height, different values of pj will be obtained starting from different peaks. An 
alternative definition of pj, which does not suffer from this problem, has been suggested in 



In the rest of this paper, we shall assume that all pockets are identical, so the prior is 
given by Eq.(27) for symmetric potentials and by Eq.(21) with Pq{X) from Eq.(28) in the 
more general case. We note that for power-law potentials like (17), inclusion of the factor 



Pj,„„,(X;j)ocp,e^(^^^)z3(X;j), 



(30) 



[28]. 
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exp(5') does not change the quahtative character of the result. For such potentials, Eq.(22) 
gives 

Z ^ exp (^'/'^x) , (31) 

where we have assumed that ^ 0*x (which is always satisfied for A <^ 1). On the other 
hand, it follows from Eq. (24) that H{4>gx,X) ~ l/4>gx, and thus 

S{X) ^ (t>lx- (32) 

Hence, the exponents in the volume factor (27) and in the ergodic factor (28) are both 
proportional to (t>'^qx-i differing only by a numerical coefficient. 

It also follows from Eqs. (24), (17) that (j)q ~ A~^/("+^). Hence, we can write 

Pp„,,(A)(xexp(CA-2/(-+2))^ (33) 

where C ~ 1 is a constant. 



VII. THE SMALL-Q CATASTROPHE 

Let us now consider what Eq. (33) implies for the prior distribution PprioriQ)- Substi- 
tuting the relation (18) between Q and A in (33), we obtain 

PprioriQ) « exp(C"g-^/("+^)) (34) 

with C" ~ 1. Assuming that A is bounded from above by A ^ 1, this equation applies for 
1. 

The distribution (34) has a very sharp peak at small Q. The exponential decay towards 
larger Q overrides any power-law prefactors, and thus the large Q catastrophe is averted. 
But now there is an even bigger problem. According to the distribution (34), small values of 
Q ~ 10~^ at the low end of the anthropic range are many orders of magnitude more probable 
than the observed value of Q ~ 10~^. Thus, we have a full-blown small-Q catastrophe. 
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This situation is not specific to power-law potentials. Both the inflationary expansion 
factor and the ergodic factor depend exponentially on the parameters specifying so 
the exponential form of PpHoriQ) is generic. 



VIII. CURVATON-TYPE MODELS 



The conclusion appears to be that the density perturbations due to fluctuations of the 
inflaton are negligible, and thus the observed perturbations must have been generated by 
some other mechanism. This conclusion has been reached in [7], where it was argued that 
the fluctuations are probably due to topological defects. However, CMB observations have 
ruled out defects as the dominant source of perturbations. 

Another mechanism, which is consistent with the data, is the curvaton model [29,30]. In 
the simplest version, it involves a light scalar field x (the curvaton) of mass m <^ H. Quan- 
tum fluctuations during inflation randomize the values of x, with a Gaussian distribution of 
width ~ H'^/m. This distribution is nearly flat, dP{x) oc dx, in the range \x\ <S H'^/m. The 
energy density of the curvaton is negligible during inflation, but at a later stage it comes to 
dominate the universe. The curvaton density perturbations are of the order 



where 5x ^ H is the quantum fluctuation of x during inflation, at the epoch when the rele- 
vant scale crossed the horizon. Eventually, the curvaton decays, and the density perturbation 
(35) is transferred to the radiation. 

An alternative class of models [31,32] assumes that x is a modulus, whose value de- 
termines the masses and/or couplings of some other particles. The resulting perturbation 
amplitude is still given by (35). In both cases, we can write 



To evaluate the prior distribution for Q, we shall consider the simplest model, where the 
inflaton potential is fixed. If, for example, the potential is given by Eq.(17), then we assume 



(35) 



Q ~ H/x. 



(36) 
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that A has the same value everywhere in the universe. The variation of Q is then entirely 
due to the long- wavelength fluctuations of the fleld x- 

The crucial difference from the case of inflaton-induced perturbations is that the energy 
of X is subdominant during inflation, and thus the inflationary expansion factor Z is not 
correlated with Q. The prior distribution for Q is 

dPprioriQ) oc oc dQ / , (37) 

and Eq.(15) gives 

dP{Q) oc QdQ. (38) 

This appears to be an improvement over the models that we considered so far. The 
exponential dependence on Q has disappeared, and the growth of the probability towards 
larger Q is considerably milder than before. However, the value of Q ~ 10~^ is still 100 
times more probable than the value we observe. We shall now discuss some ways to address 
this milder version of the large-Q catastrophe. 



IX. EVADING THE LARGE-Q CATASTROPHE 

A possible attitude is simply to accept that the observed value of Q has probabihty of 
about 1%. This is somewhat uncomfortable, but perhaps not unreasonably small. 

Alternatively, we can consider a multi-component curvaton model with the potential 

V{x) = \m\\ (39) 

where 

N 

X'^Y^XI- (40) 

a=l 

In this case, the distributions (37) and (38) are replaced by 

dPprior{Q) OC X^'-^dx oc Q-^''+^UQ (41) 
15 



and 

dP{Q) (X Q^-^dQ. (42) 

For N — 3, the distribution is flat on the logarithmic scale. 

Finally, suppose the distribution indeed rises towards larger Q. Then we expect to have 
Q ^ Qmax, where Qmax IS the upper boundary of the anthropic range, while the observed 
value is Q <^ Qmax- It is possible, however, that Q 

max ~ 10 ^ is an overestimate, and the 
actual bound is not so far from the observed value of Q ~ 10~^. 

The original argument by Tegmark and Rees [21] is that, as Q gets larger, galaxies form 
earlier, so typical galaxies are smaller and denser. As a result, the rate of stellar encounters 
that can disrupt planetary orbits gets larger, and becomes unacceptably large for Q ^ 10~^. 

However, it is conceivable that factors other than planetary orbit disruption play the 
key role in anthropic selection. One possibility is that perturbation of comets could be such 
a factor [17]. Comets move around the Sun, forming the Oort cloud of radius ~ 0.1 pc. 
Whenever the cloud is significantly perturbed by a passing star, a rain of comets is sent 
to the interior of the Solar system. Occasionally the comets hit the Earth, causing mass 
extinctions. Another threat is posed by nearby supernova explosions (see, e.g., [33] and 
references therein). Without attempting to estimate the rate of mass extinctions due to 
these effects, we note that the time it took to evolve intelligent beings after the last major 
extinction is comparable to the typical time interval between extinctions (~ 10^ years). This 
suggests that we are indeed close to the boundary of the anthropic range. A substantial 
increase in the rate of extinctions might interfere with the evolution of observers. 

X. CONCLUSIONS 

We studied the probability distribution for the cosmological constant A and the density 
contrast Q in models where both of these parameters are variable. With only mild assump- 
tions about the prior probabilities, the distribution P(A, Q) factorizes into two independent 



16 



distributions for the variables Q and y — h./Q^pm- (Here, pm is the density of nonrelativistic 
matter.) The distribution for y is rather model-independent and is in a good agreement with 
the observed value of y. Thus, despite recent claims to the contrary, the successful anthropic 
prediction for y is not destroyed by allowing Q to vary. 

The form of the probability distribution P{Q) depends on the assumed model for the 
generation of density perturbations and on the mechanism responsible for the variation of Q 
from one part of the universe to another. We first considered models where the perturbations 
are due to quantum fluctuations of the inflaton field 0, and their magnitude varies due to the 
variation of the inflaton self-coupling A. (The variation of A may be induced by fluctuations 
of some other light field. Alternatively, A may take a dense discretuum of values in different 
vacua of the lanscape.) In this case, P{Q) tends to have an exponential dependence on Q. 
The reason is that the same inflaton coupling A determines both the magnitude of Q and the 
amount of inflationary expansion. The probability P{Q) is proportional to the expansion 
factor, which depends exponentially on A, and therefore on Q. For power-law, "chaotic" 
inflation models, we found that P{Q) has an extremely sharp peak at small Q. Models of 
this kind cannot account for the observed density perturbations. 

We next considered curvaton-type models, where the density perturbations are generated 
by one or several light fields Xa, other than the inflaton. In this case, the magnitude of Q 
varies due to the fluctuations in the long- wavelength component of Xa-i so there is no need 
to assume any variation of the couplings. The crucial difference from the case of inflaton- 
induced perturbations is that the energy of Xa is subdominant during inflation, and thus 
the inflationary expansion factor is not correlated with Q. We find that, as a result, the 
distribution P{Q) has a rather mild dependence on Q and can even be flat in some models, in 
which case there is no conflict at all with observations. Both large and small-Q catastrophes 
can thus be averted. 

The models we have discussed are the simplest ones, with only two variable parameters. 
On the other hand, in the context of string theory landscape, we expect to have a multitude 
of variables, and perhaps none or very few parameters fixed. This multi-parameter space is 
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awaiting to be explored. Here, we shall only briefly comment on some issues relevant to the 
variation of Q. 

A simple extension of the parameter space is achieved by combining the two models 
that we have discussed: the inflaton potential has a variable coupling A, and there are, 
in addition, some light curvaton fields. As before, the probability distribution will grow 
exponentially towards small A. This growth will eventually be cut off, due to constraints 
related to thermalization and baryogenesis. The resulting distribution will have a sharp 
peak at a value of A <^ 10"^^, so most observers will find A to be close to that value. The 
density perturbations due to the inflaton will therefore be negligible. With A now flxed, the 
probability distribution for Q will be the same as in curvaton-type models that we discussed 
earlier. An extension to more complicated inflaton potentials, specifled by several variable 
parameters, requires further investigation. 

This work was partially supported by the grants FPA 2004-04582-C02-02 and DURSI 
2001-SGR-0061 (J.G.) and by the National Science Foundation (A.V.). 
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